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We theoretically analyze some of the anomalies of the optical sumrules in the high-temperature 
superconductors. In particular we address the particularly strong dependence on temperature of 
the sumrule in the normal state. Both electron-electron correlations and the presence of a Van- 
Hove singularity have been shown to enhance such a dependence. Here we consider both effects 
simultaneously by means of Dynamical Mean-Field Theory for a two dimensional Hubbard model 
with realistic parameters for different cuprates, and we find that the two effects are not cooperative, 
as they appear to compete one another in the region of parameters relevant for the experiments. 
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I. INTRODUCTION 

The phenomenon of high-temperature superconductiv- 
ity in strongly correlated materials such as the cuprates, 
still represents one of the most challenging topic in con- 
densed matter physics twenty years after the discovery. 
A reason for the difficulty in dealing with the cuprates 
is certainly the richness and complexity of their phe- 
nomenology, which gives rise to many competitive expla- 
nations, making it extremely hard to identify the relevant 
physical processes underlying the outstanding properties 
of these compounds. In this light, very important infor- 
mation would be inferred extracting from the whole body 
of experimental observations, the data measuring general 
and fundamental properties (e.g., thermodynamic quan- 
tities and sum rules), focusing as much as possible on 
ubiquitous (material independent) aspects. 

One important example, in the context of the infrared 
spectroscopy experiments, is represented by the analy- 
sis of the optical sumrule (SR) 1 . A big effort has been 
devoted in the last few years to the evaluation of the fre- 
quency integrals of the optical conductivity <j(iS) in sev- 
eral cuprates (e.g., Bi 2 Sr 2 CaCu 2 8+x (BSCCO)^^, 
La 2 - x Sr x Cu0 4 (LSCO)£, YBa 2 Cu 3 6+x (YBCO)M, 
and more recently also i/g-based cuprates^), and to the 
study of their behavior as a function of temperature and 
doping. One of the most general results emerging from 
these studies is the strong temperature dependence of 
the optical spectral weight (or partial optical sum rule) 
defined as 



temperature T c , where even the sign of the temperature 
variation of W op t depends on doping^, the main results 
of the infrared estimate of the partial SR in the normal 
phase (T > T c ) of the cuprates are essentially two: (i) A 
strong enhancement of the W op t {T, x) is observed when 
T is reduced. Such enhancement is generally quadratic in 
temperature and is not strongly doping dependent; (ii) 
The extrapolated spectral weight at T = 0, W op t{0,x) 
displays, instead, a clear doping dependence, increasing 
monotonically with x. The qualitative results do not de- 
pend significantly on the value of the cut-off £lc an d they 
hold in particular when it reaches the plasma frequency 
(around 1 eV). 

In this work we present a detailed analysis of the be- 
havior of the optical sumrule based on Dynamical Mean- 
Field Theory, a non perturbative many-body approach 
which allows for an accurate treatment of strong corre- 
lation effects, which turn out to be crucial to account 
for the experimental observations. We base our analysis 
on the previous work of Ref. ITol which we extend by 
including realistic two-dimensional bandstructures with 
Van Hove singularities. 

The paper is organized as follows: In Sec. II we dis- 
cuss previous results on non interacting models, and on 
the qualitative effect of correlations. Sec. Ill contains a 
brief introduction to DMFT and our solution. Sec. IV is 
devoted to the discussion of our DMFT results, while a 
simple interpretation of the same results is presented in 
Sec. V. Sec. VI is dedicated to concluding remarks. 
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where £lc is an upper cut-off, whose role will be discussed 
in the following. 

More in detail, leaving aside the behavior of the su- 
perconducting phase which establishes below the critical 



II. NON INTERACTING MODELS 

Most of the interest in the optical SR comes from its 
relation with the kinetic energ y 11 ! 12 . The most straight- 
forward identification is actually limited to tight-binding 
(TB) models with nearest-neighbor hopping, where, both 
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in the presence and in the absence of electron-electron in- 
teractions we have 
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where a is the lattice constant, V = a 2 the cell volume, 
and Ekin the kinetic energy of the carriers. The relation 
holds if the cutoff Sic is chosen large enough to contain 
the whole optical spectrum. For a non interacting system 
this would correspond to integrate up to the "plasma fre- 
quency" of the model. We will see briefly that this iden- 
tification can be significantly modified by the inclusion 
of next-nearest (and further) neighbor hopping. 

Eq. ([2]) , connecting the SR to the kinetic energy, leads 
to a potential trivial explanation of point (i). The T 2 
behavior can indeed be recovered through a simple Som- 
merfeld expansion of which gives indeed W^ p f = 
W? p f (0) - BT 2 , with W£B (0) cx t and B cx t^ 1 (for a 
fiat DOS of bandwidth 8i, e.g., B = e 2 7r 3 /48i, see Ref. 
0) 



On the other hand, as pointed out in Refs. lril flolfl3L this 
one-parameter model can not account simultaneously for 
the experimentally measured values of W op t(T = 0) and 
of B. In particular, realistic hopping parameters are not 
compatible with the size of the observed temperature 
variation of W op t- The experimental value of B could 
be recovered only invoking a value of t ~ 20 meV smaller 
by more than one order of magnitude than the values de- 
termined either theoretically by means of bandstructure 
calculations^, or experimentally through photoemission 
measurements^. The simple nearest-neighbor non inter- 
acting model appears even poorer when considering the 
behavior of the zero temperature sumrule (ii). First, the 
above mentioned value t ~ 20 meV is totally incompat- 
ible with the experimental values of W op t(0, x) (which 
ranges from 200 4- 500 meV). Second, the simple model 
predicts W opt (0, x) decreasing with the doping level, just 
opposite to the observations. 

A first natural step to heal the inadequacy of the 
above nearest-neighbor model is to include a more realis- 
tic bandstructure. In particular as we already mentioned, 
Eq. |2]) is valid only for pure nearest-neighbor hopping. 
When releasing this restriction, the generalization of Eq. 
© reads 
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where the sum is performed over all the Nk momenta k 
of the first Brillouin zone, ek is the dispersion and n CT (ek) 
the occupation number for a given k state. 

It has been pointed out in Ref. [l6| that in the non- 
interacting case the inclusion of the next-nearest neigh- 
bor hopping term t can determine remarkable changes in 
the above picture. More specifically, for realistic values of 
the hopping parameters, remarkable differences between 
the kinetic energy and the optical spectral weight be- 
havior (defined in Eq. ([3])) appear, mainly for dopings 



close to the two-dimensional Van Hove singularity (VHS): 
while the temperature dependence of E^in is only weakly 
affected by the doping level, W„ p f is more sensitive to 
the VHS, whose proximity determines a stronger tem- 
perature dependence. Although a stronger T-dependence 
of the spectral weight goes certainly in the direction of 
the experimental evidence, many inconsistencies remain: 
First of all the enhancement of the T-dependence of W^f f 
is still not enough to account for the experiments; sec- 
ondly for some (and realistic) values of the hopping pa- 
rameters the results can even show a change of sign in the 
temperature variation of Wj p ^ , which has never been ob- 
served experimentally. Finally -as already noted in Ref. 
Kdt - the doping dependence of WE,f(Q, x) remains oppo- 
site to the experimental data. 

The failure of these simple noninteracting models to 
capture the experimental behavior can not be surprising 
in light of the unquestioned role of electronic correlations 
in the cuprates. In Ref. we have shown that the inclu- 
sion of correlations determines indeed a huge step ahead 
in the understanding of experimental data. A strong 
suggestion to proceed in this direction comes also from 
point (ii) (see Introduction), because the monotonically 
increasing values of W op t(0,x) would find a very natu- 
ral explanation in a strongly correlated scenario, where 
the electronic mobility is minimal at half-filling and it 
increases with doping. In Ref. [l^ it is demonstrated how 
the presence of strong interactions can actually determine 
a separation of the energy scales controlling W op t(0,x) 
and its temperature dependence. More specifically, if 
one performs the frequency integral of o~(lu) up to a cut- 
off which includes both the Drude and the Mid Infrared 
(MIR) contribution, mimicking the experimental situa- 
tion, a strong T-dependence of W opt , and a qualitatively 
correct behavior of W op t{0,x) are obtained. 

On the other hand, some discrepancy with the exper- 
imental observations is present also in the data of Ref. 

m 

the DMFT data, which are computed for a Hubbard 
model with a semicircular dispersion which has no VHS, 
predict a stronger temperature dependence of W op t (T, x) 
at small doping than for the overdoped compounds, as 
a result of a larger distance from the Mott transition, in 
contrast with experiments, where this effect is not seen, 
and even an opposite behavior occurs in BSCCO^£. One 
of the reasons for such a discrepancy is that DMFT ne- 
glects spatial correlations. The effect of this neglect is 
expected to be stronger for smaller dopings 17 . A conse- 
quence of the lack of spatial correlation is the vanishing 
quasiparticle renormalization factor Z (renormalization 
of the coherent electronic bandwidth) when x — ► in 
contrast with photoemission data, in which Z is finite 
for any doping. 

Thus, a crucial step for a proper analysis of the small 
doping region is to consider cluster extensions of DMFT 
able to capture at least short-range correlations (such 
as Cellular DMFTi£ and DCAi2), where the Mott tran- 
sition can take place with a finite Z. Present Cluster- 
DMFT studies have been mainly dedicated to the su- 
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perconducting phase, which is not addressed in this 

2D 21 

pape r ' . 

In this work we do not address the cluster extensions of 
DMFT, while we try to supplement our previous single- 
site DMFT analysis by including the effects of a more re- 
alistic two-dimensional bandstructure displaying a VHS, 
which, as we commented above, can introduce remark- 
able effects for non interacting systems. Keeping in mind 
the limitations of single-site DMFT, we will not consider 
extremely small dopings, where non local correlations will 
become crucial. 

We notice in passing that the T 2 behavior of the SR 
is not limited to noninteracting systems, but it is rather 
characteristic of all Fermi- liquid systems, either weakly 
or strongly correlated. The actual coefficient is gener- 
ally a function of Z and of the imaginary part of the 
self-energy (proportional to T 2 in a Fermi- liquid) . What 
correlation effects and interaction with bosonic modes 
can change are on one hand the actual values of the 
T = value and of the coefficient B, as we have already 
found, and on the other hand they can introduce coher- 
ence scales which set the limit under which the system 
behaves like a Fermi liquid. Beyond that scale the tem- 
perature behavior can appear different from quadratic, 
but such a deviation is not necessarily the consequence 
of deep changes in the groundstate properties. 

III. DMFT FOR THE TWO-DIMENSIONAL 
HUBBARD MODEL 

In this paper we consider the two-dimensional Hubbard 
model, i.e., 

H = X!',,c>rT + (mi - \) (nn - \) + 

ija i ^ ^ ^ ^ 

i 

where Ci CT (c| (T ) are annihilation(creation) operators for 
fermions of spin a on site i, ni a = cj^c^, and the sum in 
the first term includes nearest-neighbors (NN) and next- 
to-nearest neighbor (NNN) hopping processes, whose am- 
plitudes are given by —t and t' respectively . 

Our choice of the parameter appearing in the Hamil- 
tonian ([!]), namely t,t',U aims to the closest contact 
with the cuprate properties. Specifically, we concentrate 
on two of the most studied cuprates, i.e., BSCCO and 
LSCO, for which a certain agreement about the estimates 
of the hopping parameters has been reached both on the 
theoretical (density-functional theory calculations) and 
the experimental side (photoemission) . In particular, for 
both BSCCO and LSCO the NN hopping parameter t 
is estimated around AOOmeV, while the different crys- 
tal structure of their unit cells reflects in different NNN 
hopping term: LDA calculations predicts a t' = 0.17t 
for LSCO, whereas larger values of t' are estimated for 
BSCCO t' = (0.25 0.30)t. Let us just note that, in 



the non-interacting case, these values correspond a VHS 
located at lower energy with respect to the half-filling 
chemical potential: In the case of t' — 0.17t (LSCO), the 
VHS would be crossed at a doping level slightly below 
x = 0.15, while for = 0.30t (BSCCO) the VHS would 
not be reached even for the higher doping levels relevant 
for the cuprates, being located roughly at x = 0.29. 

The choice of the value of the repulsive term U is cer- 
tainly less obvious: Although the relevance of electron- 
electron correlations is widely, if not universally, accepted 
as a key element in the properties of the cuprates, it is 
difficult to estimate precisely its value. In this paper, we 
have chosen U = 12t, consistently with Ref. [l(| which al- 
lows for a direct comparison with the results for the sim- 
pler semicircular DOS used in that paper, and it allows 
for a clear separation between the MIR and the Hubbard 
contributions in the optical sum rule (see next section). 
More generally, the choice of U = I2t is also guided by 
experimental evidence from the neutron scattering data 
in the cuprates, which estimate the antiferromagnetic ex- 
change J = At 2 /U ~ 140 meV. A recent DMFT analysis 
suggests however that a slightly smaller value U = lOi is 
more appropriate to describe the zero temperature opti- 
cal spectra and their doping dependence. 22 . 

As we anticipated, we use DMFT to solve the Hamilto- 
nian (|3} . DMFT is a nonperturbative method that maps 
a lattice model onto an effective local model, in which the 
effect of the neighboring sites on a given site is expressed 
through a "dynamical Weiss field". The local effective 
model can be parameterized by an Anderson Impurity 
Model (AIM), in which an interacting site is hybridized 
with a non interacting bath descri bing the Weiss field. 
For more details we refer to Ref. [23l The mean-field 
is enforced by a self-consistency condition which relates 
the Green's function of the effective model to the local 
component of the lattice Green's function of the original 
model. Namely 

where G(ui) and S(w) are the Green's function and the 
self-energy of the effective model and £k is the bare dis- 
persion on the chosen lattice (from now on, we are set- 
ting a, e = 1) . The self-consistency condition requires 
to solve iteratively the AIM until Eq. ([5]) is obeyed. 
In this work we consider the two-dimensional dispersion 
£k = —2t(cosk x + cos k y ) + At' (cos k x cos k y ) , while in 
Ref. [13 we considered a simple semicircular density of 
states characteristic of an infinite coordination Bethe lat- 
tice. We emphasize that the lattice structure enters the 
DMFT only in the self-consistency condition. 

We use exact diagonalization (ED) to solve the AIM. In 
this method the model can be solved at zero temperature 
by discretizing the bath function into a small number of 
levels, which here will be Nf, = 7. The main limitation 
of the approach is that the spectral properties are those 
of a finite system, hence the fine details can not be re- 
solved with great accuracy. This choice is particularly 
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useful for the subject of this work, considering that we 
are interested in integrated optical spectra, rather than 
in their details, and that the temperatures of interest are 
very small with respect to the energy scales of the Hamil- 
tonian ([4J (the range of temperature < T < 30QK 
corresponds to < T < 0.065i). 

In particular, the very low temperature range allows us 
to exploit the Lanczos algorithm at finite T, which has 
been developed for DMFT in Ref. [24j. In this scheme 
we avoid to compute the whole spectrum of the Hamil- 
tonian, limiting to the relevant low-lying states, allowing 
for a faster calculation than the full diagonalization of 
the matrix. 



IV. DMFT RESULTS: SPECTRA AND 
OPTICAL CONDUCTIVITY 

Ad discussed in Ref. [13 and in several precedent works, 
the interpretation of the optical spectra in DMFT can 
be greatly helped by an inspection of the single-particle 
spectral function. In particular the k— integrated spectral 
function, namely the interacting density of states, can be 
easily computed through the knowledge of the retarded 
self-energy E re t(w) on the real axis 
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FIG. 1: Interacting density of states N(u>) computed via 
DMFT at low temperature (T = 0.02£) for the case of 
BSCCO (t' = 0.3*) and LSCO (t' = 0.17t) for two differ- 
ent values of doping. Note the presence of a narrow QP at 
the Fermi level, very close to the lower Hubbard band, while 
the upper Hubbard band is clearly visible at uj — U = 12t. 



N(u) = — V Im — — (6) 

TiN k ^ w + ^-e k -E ret (w) 

In Fig. [T] we plot N(uo) for T = 0.02i, the smallest 
temperature we considered for two doping levels and for 
parameters corresponding to the two different materials. 
Quite generally, the evolution of the density of states is 
analogous to the Bethe lattice case and contains three 
features: (i) a strongly renormalized Quasiparticle peak 
(QP) at the Fermi level (of width W = ZW, where 
Z = (1 — dY,(uj) / duj)^ 1 is the quasiparticle residue), 
which is basically attached to (ii) the lower Hubbard 
band (of width roughly equal to the bare bandwidth 8t) , 
and well separated from the upper Hubbard band, whose 
center is located at u> ~ U = 12t, again with width close 
to 8t. The weight of the QP clearly increases with the 
doping level, while the spectral gap between the QP and 
the Hubbard bands appears rather stable in the range of 
doping considered (being of order of 6t 4- 8t). The effect 
of t' is hardly visible in N(ui) and seems to affect mainly 
the shape of the QP, as can be seen by comparing the 
results for the two values of t' corresponding to BSCCO 
and LSCO. 

Information about the spectral function are particu- 
larly useful to understand optical spectra due to the sim- 
plifications introduced by DMFT in this regard. The 
locality of X(w) and of the two-particle irreducible ver- 
tices in DMFT, together with the odd symmetry of the 
current operator for k— > — k determine in fact the van- 
ishing of all the vertex corrections to the current-current 
paramagnetic kernel 23 . As a consequence cr{uS) can be 



computed via the simple "bubble" 

= ^VfcX!^ J dv ImG ret (e k ,u) 

k 

ImG ret (^,uj + v) (7) 

where v£ = (ffr;) 2 , G ret (e k ,uj) = (uj + (i - e k - 
S ret (o;)) _1 is the retarded Green function, and finally 
f(uj) = (e"/ T + 1) _1 the Fermi function. For sake of 
simplicity, the normalization of Eq. [7] has been cho- 
sen so that W op t is approaching directly the value of 

Examples of the DMFT results for a(u>) are shown in 
Fig. [5] for the same set of parameters of the DOS of Fig. 
Q] One can immediately recognize how the main features 
of N(cu) discussed above reflect in a(u>). Similarly to 
the results of Ref. 10, the optical conductivity displays 
(i) a clear Drude peak at low frequencies (uj < ZW), 
determined by optical transitions occurring at energies 
within the QP width, then (ii) a Mid Infrared (MIR) 
bump at uj ~ W/2, related to transitions between the QP 
and the lower Hubbard band, and finally a high energy 
contribution at uj ~ U, which is related to transitions 
involving the upper Hubbard band. 

The rather neat separation between the low energy fea- 
tures (Drude and MIR) and the higher energy Hubbard 
contribution suggests a value of Q c ~ 6t 4- 8t as a natu- 
ral cut-off to compare (see next section) our calculations 
with the experimental data. This choice has the same 
spirit of choosing the plasma frequency as a cut-off in 



5 



0.05 



Dnide 



BSCCO, x=0.l2 



J Hubbard 
\ 




BSCCO, i=0.26 



0.D5 




0.05 



LSCO, x=0.26 



0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



0.5 1 1.5 2 2.5 3 

co/4t 



400 
380 
360 
340 
320 
300 
280 
260 
240 



220 







¥'"' 




•"'© 


BSCCO: SR(T=0) 
LSCO:W opl (T=0) . 
LSCO: SR (T=0) 



0.05 



0.1 



0.15 



0.2 



0.25 



0.3 



FIG. 2: Optical conductivity <t(u) computed via DMFT at 
low temperature (T = 0.02t) for the same parameters of Fig. 
[T] One can immediately notice the narrow Drude peak at 
low-frequency, a broader MIR contribution at ui ~ 4t and a 
small Hubbard term at higher frequencies. 



FIG. 3: Optical spectral weight including Drude and MIR 
contributions, and total SR computed in DMFT at T — as 
a function of doping for the case of BSCCO and LSCO. 



the experiments in order to select only the contribution 
from the lowest absorption band. Yet, for our choice of 
parameters, the cut-off we use is 1.5-^2 times larger than 
the plasma frequency of BSCCO and LSCO, but using 
the experimental values in our calculations would result 
in including only a part of the MIR in the integral, mak- 
ing the analysis less significant, at least in our view. 



V. DMFT RESULTS: OPTICAL SUM RULES 

In this section we analyze the results for the opti- 
cal spectral weight W op t of the Hubbard model through 
direct numerical integration of the optical conductiv- 
ity using a cut-off 8t in Eq. ([1]). From a theoretical 
point of view, it is also interesting to compare W op t 
with the total SR (including in this case also the Hub- 
bard contributions), given by Eq. ([3]), where n r7 (e] s ) — 
— 1/tt J duj ImG(ek, u>). 

Even if our focus is the dependence on temperature of 
these two quantities, a brief analysis of the T — behav- 
ior is necessary before moving to finite temperature. In 
Fig. [3]we report the values of the extrapolated T = val- 
ues for W op t and SR for three different doping levels for 
the set of hopping parameters which refer to BSCCO and 
LSCO. Both quantities increase as a function of doping, 
analogously to the Bethe lattice case. In this respect, 
as one could have expected, considering a more realis- 
tic bandstructure dispersion reflects at most in minor 
corrections of the DMFT results of Ref [l(J, and in mi- 
nor differences between the two compounds. This clearly 
shows that the strong correlation effects play the domi- 



nant role in determining the values of both W opt (T = 0) 
and SR(T = 0), regardless the shape of the non inter- 
acting bands. We notice that, apart from the natural 
difference (which ranges from 80 -j- 40 meV) due to the 
contribution of the Hubbard bands, the doping behavior 
of W opt {T = 0) roughly tracks that of SR(T = 0), since 
for the doping considered the magnitude of the frequency 
integral in the SR is mostly determined by the MIR and 
the Drude contributions. 

We turn now to the temperature dependence of the 
spectral weight. The results of the DMFT calculations 
are reported in Fig. 2J where the ratio between the value 
of W opt {T) (and the SR) and its corresponding T = 
extrapolation for the case of BSCCO and LSCO at x = 
0.12,0.19,0.26 are shown as a function of temperature. 

We see immediately that for both BSCCO and LSCO 
the relative variation W opt (T) /W opt (0) between T=0 and 
room temperature (in our units T = 300K corresponds to 
0.064£) is larger for the underdoped compounds and con- 
stantly decreasing with doping. For the case of BSCCO 
such trend is more pronounced than in the Bethe lattice 
case^ and it is only partly due to the doping depen- 
dence of W op t(0): A more quantitative analysis through 
a quadratic fit (W(T) = W - BT 2 , solid lines in Fig. 
2]) of the data at low T clearly demonstrates that the 
main contribution to this results stems from a remark- 
able reduction of the coefficient B, which controls the 
low-T behavior. For instance, in the case of BSCCO 
B ~ 21eV _1 at x = 0.12, and it decreases down to 3eV _1 
at x = 0.26. In the case of LSCO, the temperature varia- 
tion at x = 0.19 is very close to that at x = 0.26 because 
B has a minimum for doping slightly below the VHS. 

A second important observation is in order about 
the temperature dependence of the sum rule. With 
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FIG. 4: Relative variation with respect to T = of W op t (left 
column) and of the total SR (right column) as a function of 
T for BSCCO (first row) and LSCO (second row). The solid 
lines are the results of the quadratic fit at low temperatures 
(see text). 



the only exception of overdoped LSCO, in the param- 
eter region considered here, the temperature variation 
of W op t{T) /W pt(T = 0) is always bigger than that of 
SR{T)/SR{T = 0) (for BSCCO by more than a factor 
2): This results is therefore not determined simply by the 
small difference between W opt (T = 0) and SR(T = 0): 
The coefficient of the quadratic fit is always lower in the 
case of the SR indicating that the temperature depen- 
dence of the high-energy Hubbard contribution has al- 
ways the opposite sign with respect to that of W op t- This 
should be taken into account when comparing theoretical 
prediction on the spectral weight computed in presence 
of strong interactions with the experimental data. This 
problem can be faced, e.g., in Cluster DMFT, where the 
difficulty to evaluate the vertex corrections may suggest 
to resort to the calculation of the total SR. 

The main outcome of these DMFT results is that in 
presence of a strong correlation the effects of the details 
of more realistic bandstructure are not very evident in the 
behavior of the integrated optical spectral weight. As one 
clearly sees in Fig. [3] and [4] the behavior of W op t(T = 0) 
and of its dependence on temperature when the doping 
level is enhanced from x = 0.12 to x = 0.26 is dominated 
by the quasiparticle renormalization factor Z which ob- 
viously increases with the distance from the Mott transi- 
tion. In this respect, our DMFT results show that most 
of the peculiar trends determined by the two-dimensional 
density of states in a simple non-interacting scheme^ are 
washed out by strong interaction effects, so that the main 
outcomes of Ref. USl remain unaltered. 

Nonetheless we observe some effects of the more re- 
alistic bandstructure also in the presence of strong cor- 



relation, and the outcome is not completely obvious. In 
particular, one could have expected that, when the chem- 
ical potential moves close to a VHS, the largest temper- 
ature variation of W op t would have be achieved, because 
the proximity to the VHS and the strong interaction 
could have cooperated summed up somehow their effects 
in making more pronounced the T-dependence of W op t- 
This is apparently not the case, since the VHS is reached 
at a doping level slightly below x = 0.15 for LSCO and 
larger than x = 0.26 for BSCCO, and no trace of any 
enhancement in the T-dependence of the spectral weight 
is found there. Quite remarkably, instead, one can ob- 
serve that in BSCCO the weakening of the T-dependence 
of W opt with increasing doping is more evident than in 
the Bethe lattice, where no VHS is present, case^, and 
that in LSCO the change of the temperature variation for 
doping levels below the VHS (i.e., x = 0.19 and x = 0.26) 
is smaller than expected. This clearly suggest that the 
proximity to a VHS and the strong interaction effects are 
partially competitive. It is worth underlying that no ev- 
idence of the change of sign of the temperature variation 
when the chemical potential crosses the VHS is found in 
the LSCO data for x — 0.19 and 0.26 as opposed to the 
non interacting case studied in Ref. Il6l . 

An extremely simple scheme to understand the origin 
of these partly unexpected DMFT results will be dis- 
cussed in the next section. 



VI. RESCALED TEMPERATURE 
DEPENDENCE: A STRONGLY- RENORMALIZED 
FERMI LIQUID PICTURE 

In this section we will present a simple explanation of 
the surprising interplay between VHS and strong correla- 
tions displayed in DMFT. As we discuss below, the main 
idea behind our explanation is to introduce a simple car- 
toon for the correlated system, where the SR is given by 
the value obtained for a non interacting case with bands 
renormalizcd by the factor Z obtained in DMFT. The 
starting point of our arguments is the temperature de- 
pendence of the total SR at U — 0, which was already 
analyzed in Ref. [l6|, and it is shown in Fig. [5] for the 
bandstructure parameter considered here and different 
doping levels. 

It is worth noticing that, as it was pointed out in Ref. 
HH the low temperature behavior (which can be obtained 
also directly via the Sommerfeld expansion) displays re- 
markable changes for density close to the VHS. In par- 
ticular, defining xyns as the doping corresponding to 
the chemical potential at the VHS (i.e., /x = — At'), when 
x approaches Xvhs from below one observes a strong- 
enhancement of the slope of the SR(T) at low T, while 
as soon as x > xvhs ( e -g-j % — 0.30 for BSCCO, and 
x = {). 15 for LSCO) the sign of the slope changes abrupt- 
lyto positive, in contrast to natural expectation (and 
experimental evidence), before changing again sign for 
higher doping levels. 
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FIG. 5: Relative variation with respect to T = of the SR FIG. 6: Relative variation with respect to T = of the SR, 
for a non interacting system with parameters of BSCCO and after having rescaled the T-dependence of the non-interacting 



LSCO. 



model. 



One should also note that, as a consequence of the big 
and rapid changes occurring close to the VHS, the tem- 
perature behavior of SR can deviate rather remarkably 
from the standard T 2 behavior a la Sommcrfeld. 

As we mentioned in the previous section, however, 
most of these VHS effects, clearly visible in the non- 
interacting case, are quite weakened (and, in some cases, 
they even disappear) in our DMFT results, as a con- 
sequence of the strong correlation. To understand the 
reason for this effect we observe that in the non inter- 
acting systems the anomalous behavior of the SR due to 
the VHS is always limited to a small temperature range, 
whose size becomes smaller and smaller as the VHS is 
approached. If we consider that in presence of the strong 
interaction the low-temperature physics should be mainly 
controlled by the strongly renormalized QP at the Fermi 
level (see again the DOS reported in Fig. [I]), one can rea- 
sonably expect that correlations determine a reduction 
(by a factor Z) of the size of the low-temperature region 
which is controlled by the VHS. In particular, if the QP is 
strongly renormalized, as in the case we have considered 
here (the Z evaluated in DMFT ranges between 0.10 for 
BSCCO at x = 0.12 and 0.25 for LSCO at x = 0.26), 
the effects of the VHS should be limited to temperatures 
so low to become hardly visible in the experiments and 
in the temperature range we considered in our DMFT 
results. 

We can try to test this generic idea by mimicking 
the strongly interacting Fermi liquid physics by simply 
rescaling by a factor Z the energy scale of the non inter- 
acting system. More precisely, the renormalization of the 
QP peak determines a renormalization of the coefficient 
B of the T 2 term (B — >~ B /Z because B has the di- 
mensions of the inverse of an energy) , this means that the 
scaled results can be simply obtained by replotting the 
data of Fig. rescaling the T dependence of the SR for 
a given doping level, with the square root of Z extracted 
from DMFT (i.e., T -> T/VZ). 



The results of such rescaling are shown in Fig. Even 
at first glance, one can notice that the peculiar features 
associated to the VHS are much less evident now, and 
that the rescaled data resemble in many respects the out- 
come of the full DMFT calculation. More specifically, the 
slope of SR(T) does not change sign but in a very low 
temperature regime, with the exception of the BSCCO 
at x = 0.30, which is however above the doping regime 
of our interest. The same consideration applies to the 
increasing of the slope, which is observed for x — > xy^s 
in the non interacting case, and which is now limited to a 
very tiny low-temperature regime. These results explain 
why, in the presence of a strong correlation, the effect of 
the VHS is not only small, but even opposite to what one 
can expect from a Sommcrfeld expansion at U — 0: the 
deviation from the low-T "Sommcrfeld" regime occurs 
really at very low T, partly reducing the effect of the QP 
renormalization on the overall variation of SR between 
T = and T = 300K. As a second point, one can note 
that, apart small deviations, the low temperature behav- 
ior of the sum rule can be considered to first approxi- 
mation quadratic, in agreement with DMFT calculations 
and with most of the experimental measurements. One 
should also remark that, after the rescaling, the size of 
the temperature variation of the SR becomes much closer 
(some percent going from T = to T — 300K) to that 
of the DMFT data and the experimental observations, 
while even at the doping of the VHS the overall relative 
variation of SR does not exceed 0.5%. 

Although one cannot expect this oversimplified anal- 
ysis to reproduce exactly all the DMFT results (e.g., it 
would predict a slightly smaller temperature variation 
with respect to DMFT for LSCO), it sheds light on the 
essence of the physics which determines the sum rule be- 
havior in the presence of strong correlation. More pre- 
cisely it indicates rather clearly that the main features 
of the integrated optical spectra are determined by the 
physics of a strongly renormalized Fermi liquid, and the 
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effects of the two-dimensionaf bandstructure are gener- 
ally small. 

VII. CONCLUSIONS 

In this manuscript we presented a DMFT study of 
the optical sumrules of the two-dimensional Hubbard 
model, aiming to understand the interplay of the strong 
electron-electron correlations and the properties of the 
two-dimensional density of states, namely the Van Hove 
singularity. Even if both correlations^ and the VHS 
singularity^ have been invoked to explain the anoma- 
lously large temperature dependence of the optical spec- 
tral weight, the two effects are found not to be cooper- 
ative. Our DMFT analysis shows that correlations ac- 
tually reduce the effect of the VHS on the optical sum- 
rules. The mechanism underlying this reduced effect is 
a shrinking of the temperature scale below which the ef- 
fects of the VHS are appreciable. When correlations are 
strong (namely in doped cuprates) this effective tempera- 
ture scale becomes lower than the temperatures accessed 
in the experiments. 

This means that the inclusion of a more realistic band- 
structure cannot alterate the conclusions of Ref . [l(] about 
the necessity of considering strong correlations to get the 



correct order of magnitude of observed T-dependence of 
the optical spectral weight in the cuprates (as well as 
the doping dependence of its T = extrapolation). At 
the same time, our results imply that the almost neg- 
ligible doping dependence of the temperature variation 
of W op t(T) which has been clearly observed in several 
experiments, cannot be explained within a pure DMFT 
calculation, even including the effects of a more realistic 
bandstructure. 

In this respect, non-local corrections beyond the 
DMFT level certainly play a role, since they can intro- 
duce different renormalizations for the hopping param- 
eters and the QP weight, and the physics of the Mott 
transition can become richer, including for example pseu- 
dogap features at low temperatures. The inclusion of 
such non local effects certainly represent an interesting 
challenge for future studies. 
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